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Abstract 

The isospin breaking in the nucleon isovector axial charge, g\, are calculated 
within the external field QCD sum-rule approach. The isospin violations 
arising from the difference in up and down current quark masses and in up 
and down quark condensates are included; electromagnetic effects are not 
considered. We find 5g\/g\ (0.5 - 1.0) x 10~ 2 , where 5g\ = {g\) p + {g\)n 
and g\ = [{g\) p — (g^) n ]/2. Using the Goldberger-Treiman relation, we 
also obtain an estimate of the isospin breaking in the pion-nucleon coupling 
constant, (g pp7T0 - g n mr ) / 9nn-k ~ (2 - 7) x 10~ 3 . 
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The nucleoli isovector axial charge (or the nucleoli axial vector coupling constant) is 
defined through the nucleon matrix element of the isovector axial current at zero momentum 
transfer 

(N\u W u - d W d\N) = g\ U(p) 7^75 U (p) , (1) 

where U(p) denotes the nucleon spinor. Assuming isospin symmetry, one finds (g^jp — 
— (g\)m and the value of (g\) p = 1.2573 ± 0.0028, extracted from the neutron beta decay, 
has been quoted in the literature |3J . In nature, the isospin symmetry is broken by the current 
quark mass difference as well as the electromagnetic interaction, and thus (g\) p 7^ — (fl'^)n- 
Previous studies of the nucleon isovector axial charge in the framework of external field 
QCD sum-rule method have been made by various authors 0,0]. However, to our best 
knowledge, the isospin breaking effects have been ignored in these studies. The goal of 
this Letter is to examine the difference between (g\) p and (g\) n using the external field 
QCD sum-rule approach, which has been used in studying various nucleon matrix elements 
of bilinear quark operators p HT2j| . The isospin violation is reflected in m u 7^ m<i and the 



isospin breaking in the vacuum condensates. Electromagnetic effects will not be included. 

Let us start from the correlation function of the nucleon interpolating field in the presence 
of a constant external isovector axial vector field 

U z (q) = 1 J rf 4 xe^(0|T[r /iV (x)^(0)]|0)z , (2) 



where rj^ is the nucleon interpolating field introduced in Ref. JT3 

T) P (x) = Ubc v%(x)C'y IJ ,u h (x) 757 M 4(^) , (3) 

Vn{x) = tabc [<%(x)C'y lt d b (x)] 757 / XO) , (4) 

where u a {x) and d c (x) stand for the up and down quark fields, a, b and c are the color indices, 
and C = —C T is the charge conjugation operator. The subscript Z in Eq. @j denotes that 
we are evaluating the correlation function in the presence of the external isovector axial 
vector field Z^ the correlator Eq. (0) should be calculated with an additional term 



uj»j 5 u - dj**j b d\ , (5) 



added to the usual QCD Lagrangian. The up and down quark fields then satisfy the modified 
equations of motion: 

{ip - m u - Z^)u(x) = , (6) 

{ip -m d + Z^ 5 )d(x) = . (7) 

To first order in the external field, the correlation function can be written as 

U z {q) = U {q) + Z x U x {q) , (8) 

where Ho(q) is the correlation function in the absence of the external field which gives rise 
to the usual mass sum rules. Here we are interested in the linear response to the external 
field given by U x (q). The QCD sum rules for IT A (g) differ from those for Ho(q). The 
phenomenological representation for H x {q) at the hadron level contains a double pole at the 
nucleon mass whose residue contains the matrix element of interest. In addition there are 
single pole terms which arise from the transition matrix element between the ground state 
nucleon and excited states. These later contributions are not exponentially damped after 
Borel transformation relative to the double pole term and should be retained in a consistent 
analysis of the sum rules. On the theoretical side of the sum rules expressed in terms of 
an operator product expansion (OPE) the external field contributes in two different ways- 
by directly coupling to the quark fields in the nucleon current and by polarizing the QCD 
vacuum. 

The linear response of the correlation function, IT A (g), has three distinct invariant struc- 
tures 1,| : 

II A = n 1 (g 2 )g A ^ 7 5 + n 2 (g 2 ) 7 V + U 3 (q 2 )tq p a x P 7 5 . (9) 

So, one may derive three QCD sum rules from the three invariant functions, Il 1 (g 2 ), n 2 (g 2 ), 
and n 3 (g 2 ), respectively. In principle, the predictions based on these sum rules should be the 



same. In practice, however, one has to truncate the OPE and use a simple phenomenological 
ansatz for the spectral density; thus these sum rules usually have different merits. In par- 
ticular, some sum rules works better than the others. This pattern has been seen in various 
external field sum rules P-[T2|. As discussed extensively in Ref. ||, this may be attributed 



to the different asymptotic behavior of various sum rules. The phenomenological side of the 
external field sum rules contains single pole terms arising from the transition between the 
ground state and the excited states, whose contribution is not suppressed relative to the 
double pole term and thus contaminates the double pole contribution. The degree of this 
contamination may vary from one sum rule to another. The sum rule with smaller single 
pole contribution works better. We refer the reader to Refs. PJ3|JT0HT2| for more discussions 



on this point. As pointed out in Refs. 0,0], the sum rule obtained from Ili(g 2 ) is the most 
stable one for the problem under study. As such, we shall focus on this stable sum rule and 
disregard the sum rules based on n 2 (g 2 ) and n 3 (g 2 ). 

It is straightforward to obtain the external field sum rules following the techniques given 
in the literature. To include the isospin violation effects, we retain the terms linear in current 
quark masses and isospin breaking in the condensates. Here we truncate the OPE at the 
same level as in the previous studies. The OPE result for Ili(g 2 ) in the proton case is given 
by 

n i <? = - ~T7- 4 1 M-9 - TF^ — G ~2 + "2 + m « ^ u ~2 
1D7T 4 1671^ 7i q z Stt z q z 6tx z q z 

20 , . 2 1 2 ,_ . 1 2 _j ,. 1 
9 g 4 9 g 4 6 g 4 

where (O) = (0|O|0), and \ an d K denote the linear response of condensates to the external 

field 

(0|u7 M 7 B u- <*7 M 75<i|0)z = 2Z M x , (11) 
(0\uG^Yu - dG^Yd\0) z = 2Z^k, (12) 
with G^ u = \^^ u \pG Xp . Here we have omitted all the polynomials in g 2 which vanish under 
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the Borel transformation, and neglected the responses of the corresponding isoscalar current 
to the external isovector field. The analogous result for the neutron is 



20 .-5^2 1 2 12 , . 1 

- "77" (««/0 -T + n md («"/0 X -7 + o m « * ~4 " 



i - 1 

m d (dd) — 
Q 
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g 4 9 



(13) 



The resulting QCD sum rules can be written as 

M 4 El L-^ + * L- 4 / 9 + 1 « L~ 68 / 81 + ^ f 1 - 1 - 5 4]rh aL-"* 
4 3 3 1 2 2m/ 



20 m ^ 2 r4/9 2 A 7 Sm\ , _ 1 r 

+T (i- 7) _^_2__j max _i- 
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Ml 

M 2 



-Wi/M 2 



(14) 



for the proton case and 



M 4 El L-^ + * L~ 4//9 + 1 « L~ 68 / 81 + ^ ( 1 + 2 + ^ m aL- 
4 3 3 \ 2 2ffjj 

20 / 1 \ o r4 /9 2 A 7 - ~ 1 T 

+ - ( l +7)0 » £ v.__^l + I + _j mox _ £ 
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3^ 2 2mJ X M 2 An 



-4/9 



(fid)™ , 4 



M 2 



-M 2 /M 2 



(15) 



for the neutron case, where we only keep the terms up to first order in isospin violation 
and have defined a = —(2n) 2 ((uu) + (dd) s j /2, b = (27r) 2 ((o; s /7r)G' 2 )o, m = (m u + m d )/2, 
5m = m d -m u , 7 = (dd) /(uu) -l, k = -(2tt) 2 k, x = -(2tt) 2 x, and A 2 (n) = 327r 4 A p(n ), with 
(0\i]p(n)\p{n)) = \(n) v p{n) an d ^ P (n)^p(n) = 2M p ( n ). Here A p and A n are the phenomenological 
parameters that represent the sum over the contributions from all off-diagonal transitions 
between the nucleon and excited states, and E x = \— e~ s °/ M2 + l), which accounts for 
the sum of the contributions involving excited states only, where so is an effective continuum 
threshold. We have also taken into account the anomalous dimension of the various operators 
through the factor L = ln(M 2 /AQ CD )/ln(/x 2 /AQ CD ). We take the renormalization scale \x 



and the QCD scale parameter Aqcd to be 500 MeV and 150 MeV respectively. It is 



easy to see that the sum rules give rise to (g\) p = — (g\)n when isospin violation is switched 
off. 

To analyze the above sum rules and extract the quantities of interest, we adopt the 
numerical optimization procedures used in Refs. [|IOH i"2"|j. The sum rules are sampled in the 



fiducial region of Borel M 2 , where the contributions from the high-dimensional condensates 
remain small and the continuum contribution is controllable. We choose 0.8 < M 2 < 
1.4 GeV 2 which has been identified as the fiducial region for the nucleon mass sum rules [Q. 
Here we adopt these boundaries as the maximal limits of applicability of the external field 
sum rules. The sum-rule predictions are obtained by minimizing the logarithmic measure 
5(M 2 ) = ln[maximum{LHS,RHS}/minimum{LHS,RHS}] averaged over 150 points evenly 
spaced within the fiducial region of M 2 , where LHS and RHS denote the left- and right-hand 
sides of the sum rules, respectively. Note that the coupling strengths A 2 ^ also appear in the 
external field sum rules. Here we use the experimental values for the nucleon masses and 
extract A 2 ( n ) from the nucleon mass sum rules given in Ref. Jl4| [Eqs. (17) and (21)], using 
the same optimization procedure as described above. Since we neglect the electromagnetic 
interactions, we correct the neutron and proton masses such that M n — M p = 2.06 MeV, 
where the central value for the electromagnetic contribution, (M n — M p ) c \ = —0.76 MeV |15 



has been used. We then extract (g\) P ( n ), A p ( n ), and (so) p ( n ) from the sum rules Eqs. (|l4Hl5"|). 



For vacuum condensates, we use a = 0.55 GeV 3 (m u + m d ~ 11.8MeV) and b = 0.5 GeV 4 
H. The parameters x an d k have been estimated previously. Here we just quote the values, 
X = — 2/ 2 and k ~ —0.2% NTH], with f n = 93 MeV. The quark mass difference 5m has 
been determined by Gasser and Leutwyler, 5m/ {m u + = 0.28 ± 0.03 |}I5| . The value 



of 7 has been estimated previously in various approaches []T7H^5|] . All of them, with the 



exception of Refs. p4] , |25|| , indicate an interval of —0.01 < 7 < —0.006. For these 7 values, 
we obtain 

5g\ ~ (0.5- 1.0) x 10~ 2 , (16) 
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FIG. 1. The left-hand side (solid) and right-hand side (dashed) of Eq. ( |14| ) as functions of Borel 
M 2 , with 7 = —0.008 and the optimized values for (g\) p , A p and (so)p- 

where 8g\ = (#i) P + (#i)„ and g\ = [(g\) p -(g\) n }/2. With smaller magnitude for 7, we get 
smaller values for 8g\. To see how well the sum rules work, we plot the left- and right-hand 
sides of the sum rule Eq. ( |HD as functions of M 2 in Fig. [I], with 7 = —0.008. One can see 
that the two sides have a good overlap. This is typical for other values of 7 and for the 
neutron case. 

As emphasized above as well as in the literature, the contribution from the transition 
between the ground state nucleon and the excited states is not suppressed relative to the 
double pole term of interest. This contribution is included through a single constant pa- 



rameter, A p ( n ). This, as pointed out in Ref. [pq| , is an approximation. In principle, A p ( n ) 
should also be dependent on M 2 . The impact of approximating A p ^ as a constant on the 



extracted quantities is expected to be small ||26|| . Moreover, we have treated the continuum 
threshold s p ^ as a free parameter to be extracted from the sum rules. This should partially 
account for the M 2 dependence of A p ( n y 

It is also worth pointing out that unlike the mass there are no experimental values for the 
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couplings, Xpt n y One usually evaluates these parameters from the mass sum rules by fixing 
the mass at the experimental value. This means that the uncertainties associated with Xp/ n ) 
will give rise to additional uncertainties in the determination of the nucleon matrix elements 
of various current, besides the uncertainties in the external field sum rules themselves. This 
is a general drawback of external field sum- rule approach and/or QCD sum-rule calculations 
based on three point functions. Here we have not considered the uncertainties associated 

From our calculation of the isospin breaking in the nucleon isovector axial charge, we 
may estimate the isospin splitting in the pion-nucleon coupling constants by invoking the 
Goldberger-Treiman relation. The pion nucleon couplings are defined through the interac- 
tions 

£pp-K = 9pp-K Pilb^0P , £-nnn = -Qnn-Kjtil^Qn . (17) 

Noth that both 

9ppir and g nr m are positive in this notation. The Goldberger-Treiman 
relation then states 

(9a)p ~ 9ppno i i.9A)n — ~9nmro • (18) 

Using our results for (g\)p an d (<?i)n> we hud 

9ppn - 9nn, ^ (2 _ 7) x 1Q -3 ? (lg) 
9NNtt 

where gNNwo — (9pp-k + 9nmr )/2- Since the Goldberger-Treiman relation only holds approx- 
imately and there may be corrections to this relation which are originated from the isospin 
breaking effects, our estimate here is only qualitative. Nevertheless, the estimate given by 
Eq. (|i9|) is qualitatively compatible with the recent result obtained by Henley and Meis- 



sener |27| from the QCD sum rules based on three point function, though the magnitude 



is somewhat smaller than that given in Ref. W<n - O ur estimate here is also consistent with 



those found in various models |28H3lH . On the other hand, Ref. [3^] gives a result which has 
a opposite sign. 
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In summary, we have calculated the isospin breaking in the nucleon isovector axial charge 
within the external field QCD sum-rule method. We included the isospin breaking effects due 
to the difference in current quark mass difference and in quark condensates, and neglected 
the electromagnetic effects. We found a small isospin violation in the nucleon isovector axial 
charge, 8g\/ ' g\ (0.5 — 1.0) x 10~ 2 . This, upon using the Goldberger-Treiman relation, 
leads to an estimate of the isospin breaking in the pion- nucleon coupling constant, (g ppno — 
9nnw )/gNNn ~ (2 — 7) x 10~ 3 , which is qualitatively consistent with previous studies. 

This work was supported by the Natural Sciences and Engineering Research Council of 
Canada. 
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